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ABSTRACT 

We consider combined rotational, magnetic, and thermal evolution of the neutron star during 
the accretion phase in a binary system. A rapid accretion-driven decay of the magnetic held 
decreases substantially the efficiency of angular momentum transfer. As a result, the neutron 
star cannot spin up to periods shorter than some limiting value even if accretion is very long 
and accretion rate is high. The proposed mechanism can explain a discrepancy between the 
shortest detected period and minimal possible spin period of neutron stars. 
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1 INTRODUCTION 

According to the generally acc epted point of view 
jBhattacharya & van den Heuvell Il99ih the population of ra¬ 
dio millisecond pulsars originates from low mass X-ray binary 
systems (LMXBs). The recycling scenario suggests that a neutron 
star bom in a binary system with more or less typical parameters 
can accrete the gas stripped by the companion until its magnetic 
field and spin period will fit the so-called "millisecond pulsar box” 
in the B-P (magnetic field - period) plane. The characteristic 
values of the magnetic field and period for millisecond pulsars are 
B ~ 10® — 10® G and P ~ 1 — 10 ms, respectively, and these 
values can be reached as a result of evolution in a binary system. 
The evolution of a neutron star in a low-mass binary is extremely 
long and complex, because even the main-sequence lifetime of 
a low-mass compan ion exceeds 10® yrs. I n accordance with the 
stand ard scenario llPringle & ReesI 1 19721 : llllarionov & SunvaevI 
Il975h the neutron star in a close binary passes throughout several 
evolutionary phases: 

I) The initial phase, in which the pressure of the pulsar radiation 
keeps the wind plasma of a companion away from the neutron star 
magnetosphere. The magnetic, thermal, and rotational evolution 
do not differ from those of an isolated star. 

II) The propeller phase, in which the radiation pressure reduced 
by the spin-down and field decay cannot prevent the wind plasma 
from interaction with the magnetosphere. However, rotation is still 
sufficiently fast to eject the wind plasma by a propeller effect. 

III) The wind accretion phase, in which the wind plasma falls 
down on to the surface of a neutron star and nuclear burning of 
the accreted material heats the star and, due to this, accelerates the 
field decay. 

IV) The enhanced accretion phase, which starts when the com¬ 
panion leaves the main-sequence and fills its Roche lobe. The 
Roche-lobe overflow drastically increases accretion on to the 


neutron star. During this phase, a vigorous mass transfer heats the 
neutron star interiors to a very high temperature, (1 — 3) x 10® 
K and accelerates essentially the field decay. The accretion torque 
spins up the neutron star to a short period. 

The pulsars processed in the above transformations can 
have finally the pa rameters close to those of millisecond pulsars 
jUrpin et al.lll998bl l3) even if they had initially the typical for pul¬ 
sars magnetic fields and periods. Certainly, the neutron star ex¬ 
periences the most dramatic changes in the course of the phase 
IV, and this phase determine the final state of the pulsar. During 
the enhanced accretion phase, the neutron star initially spins up 
very rapidly until it approaches the so called spin-up line corre¬ 
sponding to the accretion rate. The spin-up line in the B-P plane 
is determined by corotation at the Al fven radius. It was suggested 
by iBhattacharva & SrinivasanI (Il99 ih that, during the further evo¬ 
lution, a balance between spin-up and the rate of field decay is 
reached, so that the neutron star slides down the corresponding 
spin-up line with a rate that is determined by the field decay. 

However, the behavior of a neutron star during the enhanced 
accretion phase can be much more complicated if one takes into 
account the combined magnetic, thermal, and rotational evolution. 
Because of a high internal temperature, the magnetic field decay 
during the accretion phase can be very fast such as mass transfer is 
not able to provide a sufficient amount of the angular momentum 
to maintain a balance at the spin-up line. Therefore, a significant 
departure from the spin-up line can be expected during the accre¬ 
tion phase. In this paper, we will show that these departures occur at 
both low and high accretion rates and, owing to them, the minimum 
period exists that can limit spin-up of the accreting neutron star. The 
present study is addressed mainly the neutron stars in low-mass X- 
ray binaries (LMXBs). However, the same qualitative behavior is 
typical for spin-up of accreting neutron stars in other types of bi¬ 
nary systems as well. This concerns particularly high-mass X-ray 
binaries where the accretion rate is very high and field decay is fast 
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dUrpin et al.lll998h . In such systems, the angular momentum trans¬ 
fer can also be substantially influenced by the magnetic field decay 
and thermal evolution. The considered mechanism of spin-up can 
be the key issue in understanding the origin of various classes of 
accreting pulsars. 


2 BASIC EQUATIONS 


Consider the evolution of a neutron star during the enhanced ac¬ 
cretion phase that begins at the end of the main-sequence life of 
a companion. The magnetic, thermal, and spin evolution of the 
neutron star may be essentially affected by such accretion. We 
assume that the magnetic field is maintained by electric currents 
in the crust. Such magnetic configuration can be generated, for 
example, by t urbulent dynamo dur i ng fir st minutes of the neu¬ 
tron star life dBonanno et alj[20()^ 1200^ . The evolution of the 
crustal field is determined by the conductive properties of the 
crust and material motion throughout it. In a very strong mag¬ 
netic field, the ohmic dissipation can be accompanied also by non- 
dissipative Hall currents. These currents affect the decay of the 
magnetic field indirectly, coupling different modes and generat¬ 
ing magnetic features with a smaller length- scale than the back- 


ground magnetic field jNaito & Koiimal 

1994j). Numerical Simula- 

tions JHollerbach & Rudigeill2002L 2004 

) indicate that some accel- 


eration of the decay of a large scale field (for instance, dipole) can 
occur if the Hall parameter is very large but the effect is not signif¬ 
icant for typical pulsar fields. Besides, this effect turns to be sensi¬ 
tive to the initial magnetic geometry. For example, the presence of 
both the toroidal and dipole field components decreases the rate of 
dissipation caused by the Hall effect, and dissipation proceeds ap¬ 
proxi mately on the ohmic time scale jHollerbach & RtidigeillTOOZ 
|2004 i . Note also that the Hall currents are important only for a very 
strong magnetic fields but, at the beginning of enhanced accretion, 
the field of a neutron star is u sually weaker than the ’’standard” pul¬ 
sar field dUmin et al.lll998bl R). Therefore, we will neglect the Hall 
currents. Then, the induction equation in the crust reads 


where a is the conductivity and v is the velocity of crustal mat¬ 
ter. The velocity v is caused by the flux of the accreted matter and 
is non-zero only during the periods when the neutron star under¬ 
goes accretion. This flux carries out the magnetic field into the deep 
crustal layers. Assuming spherical symmetry of the material flow, 
the velocity in the negative radial direction can be written as 


M 


( 2 ) 


where M is the accretion rate and p is the density. We con- 
sid er the evolution of a n axisymmetric magnetic field follow¬ 
ing IWendell et alj l ll987h . Introducing the vector potential A = 
(0, 0, where = S{r,9,t)/r and {r,6,ip) are the spher¬ 
ical coordinates, we obtain the equation for S from Eq.(l). The 
function S can be separated in r and 9 in the form 


OO 

S = ^si{r,t)Pl {cos 0), (3) 

1=1 

where P/ (cos 9) is the associated Legendre polynomial with the 
index 1. Using the properties of Legendre polynomials, we have 


dsi _ dsi d^si l{l + 1) 

dt ^ dr Tttct dr^ 


(4) 


In the case of a dipole field (I = 1), the function si{r,t) can 
be related to the surface magnetic field at the pole , by 

B^jt ) = 2si {R, t)/R? where R is the stellar radius l lUrpin et al.l 
Il994h . Since Eq. (5) is linear, we can normalize si in such a 
way that si(P,0) = 1. Then, the ratio Bp{t)/Bp{0) is given by 
si{R,t). 

Continuity of the magnetic field at the stellar surface r = R 
yields the boundary condition for Eq.(5) 


dsi I _ 


(5) 


In this paper, we consider the case of a dipolar field with 1 = 1. 

Since the crustal conductivity depends on the temperature T, 
accretion must influence the field decay by changing both T and 
inward directed flux of matter. Therefore, the magnetic evolution 
of a neutron star in a binary turns out to be strongly coupled to 
its thermal evolution. The thermal evolution is determined mainly 
by pycnonuclear reactions and, due to them, the crust is heated by 
nuclear burning of the accreted material. The enhanced accretion 
heats the neutron star to a high temperature ~ (1 — 3) x 10® K that 
leads to a rapid field decay. The thermal structure of the star with 
pycnonuclear reactions has been studied by a num ber of authors 
terown & Bildstenll 19981 : iHaensel & ZdunilJl2008L and references 
therein). In LMXBs, accretion due to Roche-lobe overflow can last 
as long as 10® — 10® yrs and the magnetic field can be drastically 
reduced in the course of this phase. 

Our model is in contrast to the widely accepted assump¬ 
tion of a proportionality between to the amount of accreted ma- 
terial and the field de cay rate at any time along the spinup line 
dShibazaki et alj[l989h . This simple model is obviously inconsis¬ 
tent from a theoretical point of view because the decay is deter¬ 
mined not only by the duration of accretion (or, equivalently, by 
the total amount of accreted mass) but also by the conductivity of 
the crust. The latter is determined by the temperature which is de¬ 
pendent on the accretion rate. As a result, a consistent description 
of the field decay must take into account both the duration of ac¬ 
cretion (or total amount of accreted mass) and ac cretion rate. A de¬ 
tailed study of the accretion-driven field decay l lUrpin & Gepperd 
Il995h reveals that, generally, the decay is not proportional to the 
amount of accreted mass at any accretion rate. Departures from 
this simple relation are essential even at a low accretion rate but 
the dependence becomes particularly complicated if the accretion 
rate is high, ^ 10~^° Mp)/yr (see, e.g., Eig. 5 of the paper by 
lUrpin & Gepperd (Il995h ). It is worth noting that a simple model 
with a field decay proportional to the amount of accreted material 
is not in agreement w ith observatio nal data as well. This point is 
discussed in detail by ( IWiierslll997h . In the model developed in the 
present study, the magnetic field decay is determined by both the 
accretion rate and duration of accretion. 

The rotational evolution is entirely determined by accretion 
as well because the mass flow carries a large amount of the an¬ 
gular momentum to the neutron star. It is widely believed that the 
mass flow interacts with the magnetosphere at the so-called Alfven 
radius, Ra, which is determined by a balance of the magnetic pres¬ 
sure and the dynamical pressure of plasma. 


( 2R^Bl \ 

V4Mv/G^/ 


( 6 ) 


where M is the neutron star mass and G is the gravitational con¬ 
stant. If the neutron star rotates rapidly and its angular velocity 
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is larger than the Keplerian angular velocity at the Alfven radius, 
Q. > Q.k(Ra) = plasma penetrating to the Alfven 

radius should be provided some portion of the angular momen¬ 
tum from the rapidly rotating magnetosphere and, as a result, this 
plasma must be expelled. On the contrary, matter carrying the an¬ 
gular momentum can falls down on to the neutron star if its rotation 
is relatively slow and O < Qk{Ra)- Due to this angular momen¬ 
tum, the neutron star spins up to a shorter period. However, the 
star cannot spin up to a period shorter than the Keplerian period 
at the Alfven radius because otherwise it will work as a propeller 
and expels plasma instead of acc reting it. It has been argued by 
IShattacharva & SrinivasanI il99 ih that a balance should be reached 
in spin up and the rate of the field decay, thus the neutron star slides 
down the spin up line with the corresponding accretion rate. 

Since the interaction with the magnetosphere occurs at the 
Alfven radius, the angular momentum carried by the accreted 
plasma can be characterized by its Keplerian value at Ra, 
Ok{Ra)R\, multiplied by some "efficiency” factor 5, ^ < 1. 
This factor depends on the geometry of accretion flow. If the ac¬ 
creted matter forms an accretion disk around the neutron star then 
^ ~ 1. If accretion is close to spherical one, then ^ is likely much 
smaller. The rate of the angular momentum transfer to the neutron 
star, J, can be estimated as ~ ^M{Q,k{Ra)R\), The correspond¬ 
ing spin-up rate is 

^ ^ "S ~ (7) 

where P = [GMRp/Af’P/ixl, 1 is the moment of inertia of the 
neutron star. The star spins up until the angular velocity of the neu¬ 
tron star Q. becomes comparable to CIk{Ra)- This condition de¬ 
termines the critical period, P^q, which distinguishes the accretion 
and propeller phases. We have for P^q 


p 

J- eq 


\2ABffRfl^ 


(8) 


where B \2 = Bp/ltT^ G, M_io = M/lO“^°M0yr“^, Mi,4 = 
M/IAMq and Re = i?/10km. Eq.(8) determines the so-called 
spin-up line in the B — P plane. Likely, the accreted matter forms 
the Keplerian disk during this phase. Therefore, spin up of the neu¬ 
tron star is driven by Eq.(7) with ^ ~ 1. At the beginning of the 
accretion phase, the neutron star approaches the spin-up line cor¬ 
responding to an enhanced accretion rate and, then, it slides down 
this line. Evolution on the spin-up line can last until either the ac¬ 
cretion regime is changed or the magnetic field becomes too weak 
to maintain a balance in spin up and field decay. In the latter case, 
the neutron star must leave the spin up line and evolve in a very 
particular way. We assume that the spin period of a neutron star 
follows Eq. (7) if P ^ Psq. However, as it was discussed above, 
we have to suppose P — P^q if Eq. (7) leads to P < Peg because 
accretion must stop at such rotation and the neutron star begins to 
work as a propeller if P < Peg, As a result of accretion, we obtain 
a rapidly rotating neutron star with a low magnetic field that can 
manifest itself as a radiopulsar after accretion is exhausted. 


3 NUMERICAL RESULTS 

The thermal evolution of a neutron star is considered by making 
use of the public code developed by Page, NSCool^ in which the 


magnetothermal evolution has been taken into account by coupling 
the induction equati on with the thermal eq uation via the Joule heat¬ 
ing as discussed bv [Bonanno et alj i2014h . In particular the induc¬ 
tion equation has been solved via an implicit scheme and the elec- 
tron conductivit y has been calculated with the approach described 
[PotekhinI l ll999ll 0. Actual calculations have been performed for 
a neutron st ar of M = lAM p, based on the APR equation of 
state (EOS) jAkmal et alJll998ll . The radius and the thickness of 
the crust for this star are 11.5 km and « 1.0 km, respectively. 
For the chemical composition of t he crust, we use th e same so- 
called “accreted matter” model by iHaensel & Zdunid JtOOSIi . In 
addition we have also considered models of M — 1.4Mq and 
M — I. 6 M 0 obtained with a stiffer EOS where the barionic mat¬ 
ter is calculated using a field-theoretical models at the mean field 
level iGlendenning & Moszkowskilll99lh . In particular the exten¬ 
sion of the crust for M = I. 4 M 0 for this EOS is about 2.1 km 
with a radius of 13.8 km, at variance with the APR EOS case. 

The impurity parameter, Q, is assumed to be constant through¬ 
out the crust and equal to 0.001. The spin period at the beginning 
of the enhanced accretion phase is assumed to be Pq = 100 s in 
all runs. Note, however, that the results are not sensitive to these 
values. The magnetorotation evolution of the neutron star depends 
on the field strength and its distribution in the crust. Our models 
at the beginning of the acc retion phase are pro duced according to 
the framework discussed in lUmin et al.l l ll998ah for wind accretion 
rate 10“^“^ Mq/ji. The magnetic field at the neutron star birth is 
assumed to be confined to the outer layers of the crust with den¬ 
sities p ^ po 10 ^® g/cm® but this choice does not influence our 
main conclusions. Calculations are made for two values of the po¬ 
lar magnetic field at the neutron star birth, Pp(0) = 3 x 10^® and 
6 X 10^^ G. This magnetic field processed in the course of phases 
I, II, and III then was used as the initial magnetic configuration for 
the phase IV. Our calculations of phase IV stop when the magnetic 
field of a neutron star approaches the value Bp = 10® G because, 
at the present, it is difficult to detect pulsars with Bp < 10® G. 

In Fig. I, we plot the evolutionary tracks of a neutron star for 
three different accretion rates during the enhanced accretion phase, 
M = 2.2 ■ 10“® (solid line), 4 • 10”® (dashed), and 5 • 10"^“ 
MqIjt (dot-dashed). The critical temperature for the neutron ®P 2 
gap is 10® K and the magnetic field at the beginning of enhanced 
accretion is 8 • 10^^ G as it follows from our choice of the ini¬ 
tial condition for phase IV. After accretion starts, the neutron star 
moves rapidly to the corresponding spin-up line and after 9.8 • 10® 
yrs reaches it. In its further evolution, the star slides down the spin- 
up line until the balance is maintained between the spin-up and rate 
of the field decay. However, this balance cannot be maintained for 
a long time because the magnetic field will eventually decay and 
become too weak for that. Due to the field decay, the Alfven radius 
becomes small where the accretion flow interacts with a magneto¬ 
sphere and, therefore, the amount of angular momentum transferred 
to the neutron star is small as well. As a result, the neutron star be¬ 
gins to depart from the spin-up line. This behavior is typical for 
all considered models and can be clarified by Eq. (7). The rate of 
angular momentum transfer is proportional to Bp^ and, hence, it 
decreases with decreasing Bp. On the contrary, the rate of field de¬ 
cay does not depend on the magnetic field and, at a certain moment, 
the rate of angular momentum transfer should become smaller than 
the rate of field decay. Starting from this moment, the star cannot 
slide down the spin-up line anymore. Departures from the spin-up 
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line become significant, however, at different P for different ac¬ 
cretion rates. A further evolution leads to a formation of the pulsar 
with low magnetic field and relatively long spin period. We define 
the minimum period. Pm, that the star can reach at a given accre¬ 
tion rate as the period that a neutron star has at Bp = 10® G, and 
we stop calculations at that time. For models considered in Fig. 1, 
Pm is equal to 6.5 ms for the accretion rate of 2.2 ■ 10~® MqIji 
(solid line in Fig. 1). The minimum period is determined by differ¬ 
ent mechanisms of field decay at different accretion rates. At a low 
accretion rate, M < 2.2 x 10~® Mq/jt, a decrease of the mag¬ 
netic field is determined mainly by the ohmic dissipation (the sec¬ 
ond term on the r.h.s. of Eq, (4)). In this case, the minimum period 
Pm decreases if the accretion rate increases. At a high accretion 
rate, M > 2.2 x 10“® Mq/jx, the advective term in Eq. (4) (the 
first term on the r.h.s.) is greater than the dissipative one, and the 
accretion flow pushes the magnetic field into deep crustal layers. 
This submergence of the magnetic field leads to its rapid dissipa¬ 
tion since we assume that the neutron star core is superconductive 
and the field cannot penetrate into the core. Therefore, the narrow 
transition layer is formed near the crust-core boundary and the field 
decays rapidly in this layer because of its relatively small thickness. 
The field decay is extremely fast in this case and its rate exceeds 
the rate of angular momentum transfer even if the magnetic field 
is relatively high. An increase of the accretion rate leads to more 
and more rapid decay of the field and, therefore, departures from 
the spin-up line can manifest itself even at higher magnetic fields. 
As a result, the minimum period increases with the increasing ac¬ 
cretion rate. Hence, the accretion rate exists at which the neutron 
star slides down the spin-up line during the longest time and, dur¬ 
ing this sliding down, the neutron star reaches the minimal value of 
minimum spin periods. If the accretion rate is higher or low than 
this value, the neutron star span up by accretion will rotate slower. 
Eor the model with the APR EOS, represented in Eig. 1, the mini¬ 
mal value of Pm is ~ 6.5 ms and the corresponding accretion rate 
is « 2 X 10~® Mq/ji. Note that the results are qualitatively the 
same for the neutron star models based on other EOS but the mini¬ 
mum periods have other values. 

In Eig. 2, we plot the minimum period defined above. Pm, as 
a function of the accretion rate for the neutron star models based 
on the APR EOS and mean-field EOS. In particular the solid curve 
represents a neutron star with a critical tem perature f or the ®P 2 gap 
of 10® K, (model “a” in IPage et all i2004h see also lElgareiv et ai.l 
( Il996h ) , and the magnetic field of 8 • 10^^ G when accretion starts. 
The dot-dashed line represents the same model of the solid line 
but with Tc = 10^° K, (gap model “c”). The dashed line in¬ 
stead correspond to the same model of the dot-dashed line, but 
with a smaller magnetic field at the beginning of the phase IV, 
Bp{0) = 1.6 X 10^^. On the other hand, the long-dashed and triple- 
dot-dashed line represent the same model of the solid-line but for 
a stiffer EOS (mean-field) for 1.4 Mq and 1.6 Mq solar masses, 
respectively. In this latter case the minimum period is as low as 4 
ms. 

By comparing the dot-dashed line and the dashed line, we no¬ 
tice that the final spin period clearly depends on the initial field 
strength being higher for a smaller magnetic field, as it is expected. 
On the other hand by comparing the solid-line and the dot-dashed 
line we observe that the spin period does not depend on the su¬ 
perconducting gap at the beginning of the accretion phase if the 
accretion rate is not very high, M ^ 2.2 x 10~® MQjyi. How¬ 
ever, this is not the case if accretion is heavy and M ^ 2.2 x 10~® 
Mq /yr the resulting minimal period is longer for higher Tc. 

We can conclude that the minimal spin period is essentially 



P [s] 


Figure 1. Evolutionaiy tracks of pulsars in the B-P plane for few accre¬ 
tion rates, M = 2.2 x 10“® (solid), 4 X 10“® (dashed), and 5 ■ 10“'^® 
(dash-and-dotted) Mq/yt. The magnetic field at the beginning of enhanced 
accretion is 8 ■ G. The long-dashed represent the same model of the 
solid line but for a stiffer EOS (mean-field) for 1.4 Mq . 

determined by the physical property of the neutron star matter. In 
particular we argue that Pm depend on the EOS and can be signif- 
icantl y shorter than ~ 7 ms if an EOS stiffer than the APR one is 
used Jumin et al.ll9981^. For example, the minimal spin p eriod is 
~ 4 ms for the EOS by dGlendenning & Moszkowskjl99lh . There¬ 
fore, measurements of the minimal spin period of neutron stars pro¬ 
cessed in low-mass binaries can provide information regarding the 
EOS. Our study shows that stiff EOSs are in a better agreement 
with observational data. 


4 DISCUSSION 

We have considered a combined magnetic, thermal, and rotational 
evolution of the accreting neutron star based on different EOSs. It 
turns out that the accreting neutron star cannot spin up to a very 
short period even if the accretion rate is high and accretion lasts 
sufficiently long time. The reason for this is a strong coupling be¬ 
tween the rotational, magnetic, and thermal evolution during the 
accretion phase. At a given accretion rate, there exists the mini¬ 
mum period. Pm, that the neutron star cannot overcome. This limit 
is caused by the fact that accretion-driven decay of the magnetic 
field leads to a rapid decrease of the angular momentum transfer. 
At some stage, the rate of the angular momentum transfer becomes 
so low that the neutron star cannot maintain a balance with the rate 
of field decay in order to slide down the spin-up line. Since the field 
decay is faster, the star leaves the spin-up line and evolves into the 
region of low-magnetic and relatively long-periodic pulsars in the 
B-P plane. 

The point where a neutron star begins to depart from the spin- 
up line can be qualitatively estimated by comparing the rate of an¬ 
gular momentum transfer with the rate of field decay. The inverse 
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0.1 1.0 10.0 100.0 
Mdot [10“'° Mg/year] 

Figure 2. The minimum period, Pm that the neutron star can reach by ac¬ 
cretion as a function of the accretion rate. The solid line corresponds to a 
neutron star with APS EOS and the magnetic field at the beginning of en¬ 
hanced accretion 8 X 10^^ G and a neutron ®P 2 gap with Tc = lO^K. 

The dot-dashed line is the same model of the solid line but with = 10^® 

K. The dashed line represents the same model of the dot-dashed but for the 
initial field of 1.6 X 10^^ G. The long-dashed and triple-dot-dashed line 
represent the same model of the solid-line but for a stiffer EOS (mean-field) 
for 1.4 Mq and 1.6 Mq solar masses, respectively. 


timescale of spin-up is given by Eq. (7), P/P ~ —Pj/2 'kI. Since 
j ~ Mt}K{RA)R\ ~ M'/GMRf^ during the enhanced accre¬ 
tion phase, we obtain 


P 0,k{Ra) MR\ 

P^ 2-kI^ ^ ^ n I 


(9) 


On the other hand, the inverse timescale of the magnetic field decay 
can be estimated as B/B ~ —1/ At where At is the time from the 
beginning of the enhanced accretion (see, e.g., Urpin & Geppert 
1995). The field decay is faster than spin-up if B/B > P/P, or 


AMR\ > I{^/Q.k{Ra)), 


( 10 ) 


where AM = MAt is the accreted mass. Generally, departures 
from the spin-up line depend on the magnetic field and the mass 
and radius of a neutron star and, hence, on the equation of state of 
nuclear matter. Neutron stars with a soft EOS (with smaller R at 
given M) need to accrete a smaller mass to leave the spin-up line. 
On the contrary, stars with a stiff EOS (greater R at given M) leave 
this line only if AM is sufficiently large. 

In this paper, we define Pm as the period reached by a neutron 
star when its magnetic field becomes equal to 10® G. This defini¬ 
tion is rather arbitrary and we introduce it mainly in order to stop 
calculations when the rotational evolution becomes very slow. In 
fact, the spin period of a neutron star continues to decrease slowly 
after it reaches the value Pm- However, this decrease turns out to 
be very slow and it becomes slower during the further evolution. 
For example. Pm is equal ~ 6 ms if we stop calculations when the 
magnetic field reaches the value 10^ G in comparison with Pm ~ 7 


ms obtained for B = 10® G in Fig.2. Besides, detection of the neu¬ 
tron stars with B < 10® G is a very complicated problem for the 
present telescopes. Therefore, we stop calculations when B < 10® 
G and restrict our discussion by the properties of neutron stars that 
can be observed at the present. 

The main reason for this relatively high minimal period of ac¬ 
creting pulsars is the rapid decrease of the magnetic field which 
leads to a low value of Ra ■ Transfer of the angular momentum be¬ 
comes too slow at small Ra ■ If the accretion flow forms the Keple- 
rian disc around a neutron star, the rate of angular momentum trans¬ 
fer can be estimated as MQ,k(Ra)Ra- Note that this estimate is 
valid only if Ra R- For the considered accretion rate, this con¬ 
dition is satisfied if the magnetic field of a neutron star is greater 
than ~ 10^ G. If the magnetic field is very weak (Bp < 10^ G) 
then Ra ^ R and the rate of angular momentum transfer is of the 
order of MQ,k{R)R^■ Since oc the angular mo¬ 

mentum transport becomes very slow in this case. Despite the inner 
disc radius cannot be smaller than the radius of the inner marginally 
stable orbit, the rate of angular momentum transport turns out to be 
low even if the accretion disc is extended to this orbit. Therefore, 
the further spin evolution of a neutron star is very slow and, in or¬ 
der to reach short periods, enhanced accretion should last extremely 
long time. 

As it was noted, the minimum period. Pm, is determined by 
the accretion rate and EOS of nuclear matter. However, this mini¬ 
mum period cannot be very short even if the accretion rate is high. 
A high accretion rate leads to a very rapid decrease of the magnetic 
field. In this case, the magnetic evolution in the crust is dominated 
by advection: the accretion flow pushes the magnetic field into the 
deep layers where the field dissipates rapidly in the narrow transi¬ 
tion region near the superconductive core. A fast submergence of 
the magnetic field leads to a rapid decrease of the rate of the angu¬ 
lar momentum transfer and substantial departures from the spin-up 
line. Due to this, the neutron star does not spend much time, sliding 
down the spin-up line, and can spin up only to relatively long peri¬ 
ods. As a result, there exists the accretion rate at which Pm reaches 
the minimal value for a given EOS. The minimal value of Pm for 
the considered model with the APR EOS is « 7 ms and the neu¬ 
tron star reaches it if the accretion rate is « 2 x 10~® Mq /yr. If 
the accretion rate is higher or lower than this value, the neutron star 
can spin up only to a longer period. This limit of the spin period is 
determined by the EOS and can differ from the value ~ 7 for other 
EOS. Note that the limit is smaller for stiff EOSs and greater for 
soft EOSs. For example, the minimal period for a neutron star with 
the EOS by jGlendenning & Moszkows^ll99lh is ~ 4 ms and this 
period can be reached if the accretion rate is ~ 4 x 10“^° Mq /yr. 
Note that, for a stiffer EOS, the neutron star reaches the minimal 
period at lower accretion rate. 

The shor test spin period dete cted up to date for PSR J1748- 
244 is 1.4 ms iHessels et al.l ( l2006h . The minimal spin periods at¬ 
tainable for neutron stars and determined by balance of the grav¬ 
ity and centrifugal forces are well below this value for most of 
the proposed EOSs of a nuclear matter. A discrepancy between 
the minimal possible spin period of neutron stars and the short¬ 
est detected period can be caused by the mechanism considered 
in this paper. The calculated values of a minimal period, ~ 7 ms 
and ~ 4 ms, clearly point o ut that the EOS should be stiffer than 


that by Akmal et ^ 1 1998 |) and, li kely, even stiffer than that by 


dOlendenning & Moszkow^ll99lh . 

Our calculations show that formation of shortly periodic pul¬ 
sars (like millisecond pulsars) is possible only if the accretion rate 
belongs to some not very wide range around M ~ 2 x 10“® Mq /yr 
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for APS EOS. If the accretion rate is lower than this value, a pulsar 
in the binary system has no chance to become a millisecond one 
even if accretion lasts very long. Note that some observational in¬ 
dications (for instance, a distribution of pulsars in the B-P plane) 
also point out on the existence of the minima l accretion rate that 
can lead to a formation of millisecond pulsars jPan et 'Zll2013h . 
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